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In this paper we study the critical exponents of the quan- 
tum sine-Gordon and U(l) Thirring models in the incom- 
mensurate phase. This phase appears when the chemical 
potential h exceeds a critical value and is characterized by 
a finite density of solitons. The low-energy sector of this 
phase is critical and is described by the Gaussian model 
(Tomonaga-Luttinger liquid) with the compactification ra- 
dius dependent on the soliton density and the sine-Gordon 
model coupling constant /3. 

For a fixed value of /3, we find that the Luttinger parame- 
ter K is equal to 1/2 at the commensurate- incommensurate 
transition point and approaches the asymptotic value /3 2 /8tt 
away from it. We describe a possible phase diagram of the 
model consisting of an array of weakly coupled chains. The 
possible phases are Fermi liquid, Spin Density Wave, Spin- 
Peierls and Wigner crystal. 

PACS No: 71.10.Pm, 65.40.-f, 65.50.+m 



In this paper we study the critical exponents of the 
quantum sine-Gordon (SG) and the U(l) Thirring mod- 
els in the presence of soliton condensate. The La- 
grangian of SG model is given by 



L= j &x\^[(d t <5>) 2 -(d x <S>) 2 ] -h(3/7rd x $ 
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where /i and (3 are real parameters, mass-like and cou- 
pling constant respectively and h is chemical potential. 
SG model is closely related to the U(l) Thirring model 
(it describes its (iso)spin sector). The Lagrangian den- 
sity for the latter model is given by 
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where j° = (l/2)# aTl^ap^fi are the (iso)spin currents. 
The relationship between the two models is established 
by bosonization (this procedure and the relationship be- 
tween the operators is discussed at length in 0). The 
coupling constants gn and gj_ are related to [3 and to \i. 

The Hamiltonian of the U(l) Thirring model, which 
as a model of two species of fermions possesses central 
charge c = 2, can be written equivalently as a sum of 



two charge- and spin- Hamiltonians, each having central 
charge c = 1, i.e. H c = v c /2 J dx (: J J : + : J J :) and 

H s = ^ J dx [: JJ : + : JJ : + g l{ : J Z J Z : + (3) 

+ g±(- JxJx ■ + ■ JyJy 0] 

In the above equations, J, J and J, J are chiral charge- 
and spin-currents given as usual by J = i?^i? Q and J = 
Ra(& a f3/2)Rf3 for the right currents and similarly for the 
left ones. 

The U(l) Thirring model may be used as a continu- 
ous model of interacting lattice fermions. Here one can 
imagine two situations. First, one may use the (iso)spin 
part of the Thirring model as a model describing the 
charge sector of the Mott insulator. In the latter case 
the currents belong to the charge sector and the spin 
sector is gapless. It is natural to impose an additional 
requirement of the SU(2) symmetry in the spin sector, 
which fixes the corresponding compactification radius 
(3 S = V8TT- In the different interpretation the model de- 
scribes the continuous limit of lattice fermions with a 
gapless charge sector and an anisotropic interaction in 
the spin sector. In this case one does not expect any 
restrictions on the compactification radii. 

In both cases one has the following relationship be- 
tween the order parameter fields and the bosonic fields: 
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with /3 c Pc — 8ir and /3 S /3 S = Sir. 

To be definite, let us stick to the Mott insulator ver- 
sion. In this case the field $ of SG model (|l|) is $ c and 
the <I> S sector is gapless with (3 S = and [3 C = (3. 
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When 1 < 8ir the cosine term in Eq.(l) is relevant and 
generates a spectral gap in the low-energy sector. How- 
ever, when the chemical potential h exceeds a certain 
threshold, the system becomes gapless again [||. 

In this paper we discuss the case [3 — > 0. This case 
is seldom discussed in the literature because such limit 
cannot be achieved by short-range interactions (recall 
that for the repulsive Hubbard model (3 varies between 
\/87r (weak repulsion) and v^tt (infinite repulsion). In 
H is found that the Luttinger liquid parameter of the 
charge sector K c takes values monotonically increasing 
in the interval 1/2 < K c = ( 2 {B)/2 < 1 as it = U/t 
decreases (the asymptotic dependence of the dressed 
charge ((B) on u is also given in j|). 

One can achieve small (3 considering long-range inter- 
actions (such as unscreened Coulomb interaction which 
can be realized, for instance, in quantum wires). 

The limit of small (3 is very curious from the math- 
ematical point of view because solitons in SG model 
become in this case infinitely heavy. Thus in the ab- 
sence of solitons (that is in the insulating phase) SG 
model becomes a model of a free massive bosonic field 
with mass fi. This field can be thought about as an op- 
tical phonon. One may wonder what happens when one 
forces solitons to appear applying the chemical poten- 
tial. As we shall demonstrate, the scaling dimensions in 
the regime of small (3 are determined not by the soliton 
mass, but by the mass of the optical phonon /i. 

Let us describe our results. In the Mott insulator in- 
terpretation where the field h truly corresponds to chem- 
ical potential we get the following. The gapless (incom- 
mensurate or doped) phase of SG model has scaling di- 
mensions dependent on (3 and soliton concentration n so i 
(which is related to deviation of the electron concentra- 
tion from half-filling) . Looking at scaling dimensions of 
various operators, one can identify different regions of 
the phase diagram. 

The most singular operators are components of the 
staggered magnetization n 2 ,^, the 2fcp and the 4fcp 
components of the staggered charge density p s ^. The 
first four operators all have the same scaling dimension 
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where K c (n so \,f3) varies from 1/2 at n so i — * to /3 2 /8tt 
in the limit of large chemical potential. The Akp charge 
density component has scaling dimension 
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Therefore d n ,d4fc F < 1 in the entire area of interest 
and the corresponding susceptibilities remain singular 
throughout this area. However, starting from K c < 1/3 
(with d^k-e < d n ) the 4/c^-component is more singular. 
Its correlation function decays more slowly being also 
oscillatory with wavelength A = 2ir/Ak-p = a, the inter- 
site spacing. In this case in the presence of other chains 



the material will undergo a phase transition to a Wigner 
crystal (4/c^-Charge Density Wave). For K c > 1/3 the 
material will undergo a phase transition either of Spin 
Density Wave (SDW) or Spin Peierls (SP) type. The 
importance of the single-electron interchain tunneling is 
determined by the magnitude of the scaling dimension 
of the single fermion creation (annihilation) operator: 



(11) 



Since the operator describing single-electron interchain 
tunneling has scaling dimension 2d,p, it becomes rele- 
vant at K c > y2 — 1 . 

From the exact solution ji|,|5| the excitation spectrum 
is given by 



E = e(9) , P = 2vr / d9'a( 
Jo 
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where e(9) is obtained from the solution of the following 
Bethe-ansatz integral equation: 

e(0)+ / d0'G(6- 9')e(9') = M B cosh0- h . (13) 



M s is the kink's mass and the Fourier image of the kernel 
is 
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The soliton's mass is a function of the coupling constant 
(3 and at small values of (3 it grows like 1/r, where r here 



and in the following is the short notation for (3 
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plot of M s versus r can be found in Q . 

The value of B in Eq. ( |l3| ) is determined by the condi- 
tion e(±B) = 0, and depends on the coupling constant 
(3 and the field h. 

At the point t = 1/2 ((3 2 = 4tt) the kernel Q van- 
ishes and e(9) is given simply by the right-hand side of 
dl3|). For other values of r, Eq. ( p^[ ) can be solved nu- 
merically. A plot of e(6), for a fixed value of the coupling 
constant (3 and different values of the field h, is shown 
on Fig. 

The soliton condensate appears when the field h ex- 
ceeds the critical value, found from (13) to be h c = M s . 
The soliton contribution to the ground state energy of 
the sine-Gordon model is given by 
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The average number of solitons n and the susceptibility 
X are given by 
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FIG. 1. Plot of —e(6) on 8 and h dependence obtained 
from the numerical solution of Eq. ([L3|), only inside the in- 
terval (—B,B), is shown. The coupling constant r is taken 
to have value 0.1. 

One can also calculate the number of solitons using 
the density distribution function a(9): 



n so i = 



d6a{6) 



(17) 



with a{6) satisfying the following integral equation: 



a(0) + dd'G(0-6')a(0') =M s cosh6/2n . (18) 



At low-energies the excitation spectrum can be lin- 
earized at the Fermi points 9 = ±-B, with the Fermi 
velocity given by: 
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The boundary B is the same as the one of Eq. ( |13[ ) . 

The scaling exponents of the Luttinger liquid are ex- 
pressed in terms of a single parameter K. This param- 
eter can be calculated in two different ways. 

First of all, it is related to the dressed charge Q(B) 
which can be found from the solution of the following 
integral equation 0: 
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where G(9) and B are the same as in Eq. ([D 
parameter K is related to (,{B) by 

K = c[C(B)} 2 

up to a constant factor which, as results from ( 
(|3l|), is c = 1/2. The scaling dimensions are 
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Another way to determine K is to use the identity 
relating K to the susceptibility x an d t ne Fermi velocity 



K = ^TTXVF 



(23) 



We can use the above formula as a check of the result 
obtained by (pl|). 

First, consider the case h — > +oo where the limit 
B goes to infinity. In this limit the integral equation 
Eq. (|l3|), as 9 — > B, can be written in the following 
form 

e(9) + / d9'G(9 - 9') e(9') = M s — - h , (24) 

J-B 1 

and the derivative e'(6) = de(9)/d9, taking into account 
the boundary conditions for e(9), fulfills 



e'(9)+ / d9'G{9 -9')e'(6') = M s — 
J-b 2 



(25) 



As can be seen from Eq ( p"8| ) and Eq. (p5|), e'(9) and 
2tt(t(9) fulfill the same integral equation, as 9 — > +oo, 
meaning that the Fermi velocity in the limit h — > +oo 
goes to 1, vp — > 1. 

According to || the ground state energy has the fol- 
lowing /i-asymptotic behavior: 

S = - T h 2 /7T, n so i = 2Th/7r, X = /3 2 /(2tt) 2 , (26) 

which, for the Luttinger liquid parameter, gives 

K = /3 2 /8tt = r . (27) 

On the other hand, we can derive the same result from 
Eq. (|20|). In the limit B — > +oo this equation can be 
written as a Wiener- Hopf equation: 



C(0)+ / d9'G{9-6')((6') = 1 



(28) 



(We have shifted 9 by B, so that now we need to calcu- 
late C(0)). The solution for 



C ( ^ ) (^)= / d9e wuj ({9) 
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where 



G(+)(0)G(-)(w)(iw + 0) 



(29) 



1 + = G (+) (w)G ( - ) (w) 

with G^ (uj) being analytic in the upper (lower) half- 
plane and G(oo) = 1. The limit at 9 = is given by 



3 



i= 10 



M =h,= 146 




FIG. 2. Plot on 6 and h dependence of the dressed charge 
((d) of Eq. (po|), only inside the interval (—B,B), is shown. 
The value of the coupling constant r is kept fixed. 

C(0) = lim iwC (-) M = [gM(o)]- 1 



= [1 + G(w = 0)]- 1/2 = V^ 



Thus we get 
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which is consistent with Eq. (27). 

Near the critical line the interval (-B, B) closes to 
a point and from (|2^) the dressed charge becomes 1. 
Otherwise, for all values of h we have solved the integral 
equation for the dressed charge, Eq. (pfj[), numerically. 
A plot of ((0), only inside the interval (-B, B), is shown 
on Fig. 0. The value of ((B) decreases (increases) fast 
near B — for r < 1/2 (r > 1/2). For larger values of 
B, ((B) reaches quickly values close to its asymptotic 
value \/2t- 

In the fermionic sector corresponding to the sine- 
Gordon model the chiral components i?f,if of the In Eqs. (|34 



FIG. 3. The dependence of the Luttinger parameter K on 
the field h for h > h c . The coupling constant r is taken to 
have values 0.1 and 0.2 respectively. In the regions K < 1/3, 
1/3 < K < V2 - 1 and \/2 - 1 < K < 1/2 the system is 
in the Wigner crystal (WG), spin density (SDW) insulating 
phases and Fermi liquid phase respectively. 



with H 



= 4 ^ 



K c -2) 



(35) 



fermion operator ifj^ are given by the expressions 
R] = cxp [-iVV2(C + C -1 )0c- 



(32) 



L-\ = cxp [ iv/^2(C - C 1 )^ + (33) 
+ i^/^(C + C 1 )^ + i\/2~^ s 

where <j> c = ($ c + 6 c )/2, </> c = ($ c - 9 c )/2 and simi- 
larly for (f> s . The correlation function of the right chiral 
component is given by 



) one should not confuse r and the 
exponent 8 with the coupling constant and the argument 
of the integral equations used previously. 

From Fig. || and Eq. ( ^l|) it is clear that for a given 
value of the coupling constant, the Luttinger parameter 
K inside the incommensurate phase changes its value 
depending on the field h (or on the number of solitons) . 
Near the critical line as h — > h^(r), the Luttinger pa- 
rameter approaches 1/2 (for all the values of the cou- 
pling constant r (= /3 2 /8tt), and from Eq. @ 6 -y 1/8, 
(here and in the following we use 6 instead of C ). For 
large values of the field h, h — -> oo, K — > r and therefore 
9 — > (t — l) 2 /At. Summarizing, for K and we have 
the following limits: 
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T, (34) 



For (9 > 1 the single-fermion density of states (DOS) 
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FIG. 4. Plot of the dependence of the exponent 9 on the 
field h obtained for h > h c . As it can be seen the exponent 
6 which enters also in the single-particle density of states 
DOS, increases monotonically with increasing h, and in case 
t < t c — 0.17157 for values of h bigger than a critical value 
ho the exponent becomes larger than 1. In this region DOS 
exhibits pseudogap type behavior. In the inset is shown the 
square root dependence of 9 on h in the region close to CL. 



exhibits pseudogap-type behavior, 
happens only when r < r c , where 



From Eq. fed) this 



t c = (3-2\/2) ~ 0.17157 



(37) 

+00) we find 9 



For all values of h in the interval (h c 
numerically. 

Plots of the dependence of the Luttinger parameter 
K and the exponent 9 on the field h are represented on 
Fig. | and Fig. | . 

((3 ~ 1, /3 2 < 8tt) Behavior of K, 9, n so \, \ 
near the critical line (CL). It is interesting to have 
the h dependence of the above quantities near the 
commensurate-incommensurate transition as h — > 
This can be achieved by solving the integral equation 
( |j"3| ) by iterations. After two iterations the boundary B 
can be written as a series in powers of (h/h c — 1) 
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with coefficients which, up to the order shown, do not 
change in higher order iterations. In this way we get for 
the soliton contribution to the ground state energy: 
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and for the dressed charge: 
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FIG. 5. Plots of the dependence of the number of solitons 
on the magnetic field (chemical potential) h. The insets 
compare the numerical (dotted line) and the analytic results 
near CL. The coupling constant r is taken to have values 0.1 
and 0.2, respectively. 
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Formulae (|39|) and ( |16| ) show that (to first order) the 
number of solitons and the susceptibility have square 
root dependence and square root singularity as functions 
of h, respectively. 

Eq. © and Eq. © show, on the other hand, that 
the Luttinger parameter K decreases, from K = 1/2 at 
h = h c , as a square root function of (h/h c — 1), whereas 
the exponent 9 increases linearly having value at h = 
h c . Explicitly, K and 9 as functions of (h/h c — 1), near 
the critical line, are given by: 



A-=i-2^G(0) (A-l 
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As functions of the number of solitons 
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they (to a first approximation) are given by: 



1 G(0) 
K = - - 27r_ ^ — "-sol (n so \ — » 0) 
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(n so] -> 0) 



A plot of the dependence of 9 on the field ft,, near the 
critical line, is shown on Fig. ^. The inset, in Fig. 
confirms the square root /i-dependence of 9 in that re- 
gion. 

In Fig. | we have represented the dependence of 
the number of solitons, present in the incommensurate 
phase, on the field h for different values of the coupling 
constant P 2 /8tt. Notice in (^ ) that the constant of 
proportionality of the number of solitons on h grows for 
decreasing r. This can be observed also numerically in 
the above plots. In Fig. || we have represented the de- 
pendence of the exponent 9 on the number of solitons. 
In Fig. Q we have represented the dependence of suscep- 
tibility on the number of solitons, outside the region, 
near the critical line, where it behaves like x ~ 1/^soi- 
Both plots are obtained numerically, with h used as pa- 
rameter. 

At small r, G(0) behaves like (1/tt 2 t) ln(l/r) and the 
expansions (39) and ( (40| ) are valid only for values of h 
which fulfill 



7r 2 T/ln(l/r) > (h/h c -l) 



(44) 



For smaller values of r the expansion ([39|) holds only 
in a decreasingly small interval h above h c . For values 
of h outside this small interval a crossover to a different 
behavior takes place. In this case one has to take into 
account the presence of the breathers in the ground state 



(J3 -> 0) The limit of small r, M s t = const. At 

(3 2 < 47r in the sine-Gordon spectrum, in addition to 
solitons, breathers (kink-antikink bound states) will ap- 
pear. The mass spectrum of the breathers is described 
by M: 



2M S sin n 



TT(3 2 



Wir-2/3 2 , 

with n - integer, taking the following values: 
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When (3 — > the bound-state energy vanishes and at 
small 1 one can expand the sine in (f45|) and obtains 



M n = nnM s T 



(47) 



As (3 — > the soliton mass grows to infinity while f3 2 M s 
converges to a finite value || . The linearity in the mass 
spectrum means that the breathers instead of being in 
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FIG. 6. Plots of the dependence of the exponent 
9 = 1/8(C(B) - 2/C(B)) 2 , Eq. (§!), on the the field h and 
on the number of solitons respectively. For large h, the ex- 
ponent 8 approaches its asymptotic value, 8 — > (r — l) 2 /4r, 
shown also in the plot. The parameter r is taken to have 
value 0.2 in both plots. 

the n-th bound-state are as n unbound particles. This 
is the classical limit, where the sine- Gordon becomes a 
Klein-Gordon model. In this limit an infinite number of 
breathers will appear in the spectrum, filling the gap to 
the soliton mass M a . This problem is discussed in ||], 
and we use in the following those results. 

The boundary B, in the region r <C B <C I , is found 
to satisfy the following integral equation {1: 
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dvhi | coth( — )| [sinh 2 _B cosh 2 w — cosh 2 !? sinh 2 w] 
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which for B will give: 



B 2 ln(l/S) w 7T 
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(In the above equations (i stands for the first breather's 
mass [i = tttM s = 7rr/i c ). The contribution to the 
ground state energy of the kink condensate is: 
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FIG. 7. Plots of the dependence of the magnetic suscepti- 
bility x on the magnetic field (chemical potential) h (upper 
plot), and on the number of solitons (lower plot), near the 
critical line are shown, r is taken to have value 0.1 in both 
plots. In the lower plot \ is shown only outside the region 
where it diverges like l/n so i (see inset of Fig. ^ for r = 0.1). 

From this we obtain the equations for the number of 
solitons and the susceptibility: 
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If we express h c /(h — h c ) in terms of the number of 
solitons, 
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TTfl 
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for the susceptibility we get: 
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These results agree with those obtained by jll],(|]. The 
criterion r <C B <C 1 means therefore: 



In order to calculate the Luttinger parameter in this 
region, we would need to have in addition to x the value 
of the Fermi velocity v F . We achive that through the 
following calculations. 

The excitation spectrum in this region is given by B 



e{9) = ~p^fsmYi{B - 6) sinh(B + 9) , (55) 
whereas the equation for the momentum is: 
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After the substitutions 

x = tanhw , Y(x) — coshvp(v) , 

we get 
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P / = g(x) 

J-b x-y 

where b = tanh B and 
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The general solution of such equation is 



Y(x) = \y/b 2 -x 2 P { 
J-} 



g(y) dy 



-b (y - x)^Jb 2 - y 2 

The solution at small b w B <C 1 is 
4n so i 



Pf ~P 



-\/b 2 — x 2 , e{x) = —pyb 2 



which gives the following dispersion law: 

E = v F (p-p F ) , 

with 
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The maximal value of — e is psmhB w pb. This means 
that E{p) deviates from the linear form at p ~ p F in- 
dicating that the ultraviolet cut-off in the momentum 
space is of order of p F ~ n so b 

According to the above results in the region (|54| ) we 
have 



v 1 

K = -n X v F « - 



7™sol 

a 



3/2 



exp 



TTfJ 



An. 



sol 



When n so i approaches the boundary of validity of this 
expression given by Eq. ( |54| ) , which corresponds to very 
small n so \/fi, K is still small: K ~ (n so i/7r/i) 2 /2. At 
even smaller n so i it asymptotically approaches 1/2. In 
this region Eq. (|2|) can be solved by iterations. 

Thus at t <C 1 (the classical regime) the Luttinger 
parameter K decreases quickly from 1/2 towards its 
asymptotic value r. This rapid change takes place for 
n so i < fi, that is, it is determined not by the soliton 
mass which goes to infinity at (3 2 — ► 0, but by the mass 
of the breather which remains constant in this limit. 

DISCUSSION 

In this paper we have studied the scaling exponents 
of the isospin sector of the U(l) Thirring model in the 
incommensurate phase when its energy spectrum is gap- 
less. This model may be used as an effective theory for 
a low-energy limit of fermionic theory with long range 
repulsion on a lattice close to half-filling. Exactly at 
half-filling the charge sector is gapped and is described 
by the sine-Gordon model with a relevant cos term. This 
gap is suppressed by a finite chemical potential and the 
system goes into a gapless incommensurate phase. In 
this phase, which is the phase of concern in this paper, 
the correlation functions of the operators (f§||) and ( pi) ) 
decay as power laws of some critical exponents which 
are given in terms of a single parameter K, the Lut- 
tinger liquid parameter. K is related to the dressed 
charge which satisfies an integral equation characteris- 
tic for exactly solvable models. 

We are mainly interested in the region of small /3. 
This is an interesting limit which has not been much 
discussed in literature. At j3 — ► the solitons become 
infinitely heavy and in the insulating phase (absence 
of solitons) SG model becomes a model of free massive 
bosonic field with the mass fi. It turns out that even in 
the incommensurate phase when the chemical potential 
forces solitons to appear, K does not depend on the 
soliton mass but on the mass of the optical phonon \x 
(see Eq. ©). 

As we have shown, the value of K depends on doping 
(that is the soliton density) with K = 1/2 in the limit of 
zero doping and reaching its asymptotic value (3 2 /8w at 
large doping. Therefore changing the number of solitons 
one can explore different phases of the model (see Fig. 
3). Thus for small doping the single-particle hopping 
is relevant. In this case one might expect that the sys- 
tem of many coupled chains will become a Fermi liquid. 
When doping increases K decreases, the single particle 
tunneling processes are suppressed, but the spin suscep- 
tibility becomes singular. This is a region when inter- 
chain coupling will establish Spin Density Wave (SDW) 
ordering. For K < 1/3 the 4fcp-component of the charge 
density operator has the smallest dimension with the 
corresponding susceptibility being the most singular. 



This notifies that the Wigner crystal order parameter 
will be stabilized by a Coulomb interaction in a three 
dimensional array of coupled chains. 

The number of solitons and the susceptibility are 
found analytically and numerically and display, at the 
transition point, square-root behavior and square-root 
singularity respectively. As (3 — > 0, the above depen- 
dence holds only in a decreasingly small interval h above 
h c . For larger h the dependence of n so i and % on /i is 
given in Eq. (51). Plots of the dependence of n so i and 



X on h are given on Fig. [| and Fig. ^. 
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study this problem. A. M. T. acknowledges the hos- 
pitality of Isaac Newton Institute where part of this 
work was completed. E. P. acknowledges discussions 
with Tilo Stroh. 

APPENDIX 
Expansion near the critical line 



The kernel G(9) of ( |13[ ) can be expanded about 9 = 
in a Taylor series with convergence radius rg = ttt/(1 — 
r). For small h c — h = M s — h, using the Taylor expan- 
sion of the inhomogeneity h c cosh(#) — h, the integral 
equation can be solved for a given accuracy by a power 
series about 9 = implying also a series in B about 
B = 0: 



N 



G (0)~El%^ > G 2l = G^(0) 



i=0 



(2i) 



N 1 

h c cosh(6>) - h ~ h c V —^9 21 - h 
w ^ {2i)\ 



(63) 



(64) 



N 



2(JV-i+l) 

Beginning with e 2 ij = the iteration for the coef- 
ficients converges after 27V + 1 cycles, where actually 
N = 4 was used and results for N = 2 are shown below. 
The solution for B from the condition e(B) = 0, where 
e(B) is now a sole power series in B, is parameterized 
as a power series in y/h/h c — 1, cf. (38), 



2JV+3 



B~ oc^ , x = ^--l 

3 = 1 

2N+2 2N+4 

e(B) ~ £ ^ ~ £ , 



(66) 
(67) 



i=0 



i=0 



where also h is replaced by h c (l + x 2 ) in e\ B \ which 
yields e[ h) = e ( Q h) = and e ( 2 h) = h c {a\ - 2)/2. From 
e !+i = 0, i = 2, . . . , 2N + 3, the coefficients a, can be 
read off recursively. 
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In the soliton contribution to the ground-state energy 
( |l5| ) the power expansion of cosh# and ( |65| ) is used, 
where the integration leads to a power series in B and 
replacement (|6|) to a series in ^Jhjh c — 1: 



So 



2V2M S 2 
3tt 



2N+2 



£3 



X = \l— - 1 
/lc 



(68) 



The integral equation ( |20[) for £ is treated exactly as 
(|l3|), leading to a series in and B. For C(-B) applying 
the previous replacement for B the series reads 



2JV+2 



c(s) - £ c 



i=0 



x = \l-r- -1 



(69) 



Finally, this series is replaced into (|2l]), which then is 
inserted into (E55j), leading to the following series: 



27V+2 

#~ £ ^ 

i=0 



2AT+2 

£< 



(70) 



£^ 6) = 1/384 - 2321G§/4200 + 286Go=/27 + 
1478048G(j/729 + 16(61/175 + 
2288Gg/27)GoG 2 /3 + 1232G^/225 + 352GoG 4 /105 



Co 

c[ h 

Jh 

Jh 
^3 

Ah 

A h 



Jh 
^6 



Coefficients C- h) of C( B ) in © 

= 1 

= -2\/2Go 
= 32G^/3 

= V2(Go/2 - 280Gg/3 - 8G 2 )/3 

= 32(-Gg/5 + 160G^/9 + 14G G 2 /5)/3 

= V2(-3G /80 + 154Gg/15 - 16016G[j/27 
2(1/3 - 336Gg/5)G 2 - 16G 4 /15) 

= 16(13G§/35 - 448G^/9 + 57344Gg/27 + 
8(-29G /35 + 704Gg/9)G 2 + 58Gi/5 + 
272G G 4 /35)/9 



Coefficients K t of K in (70) 



Coefficients on of B in (66) 



a\ = \/2, 

a 2 = -4G /3, 

a 3 = V5(-l/12 + 20Gg/9), 

a 4 = 2G /5 - 256Gg/27 - 8G 2 /5, 

a 5 = V2(3/160-49G§/45 + 616Gg/27 + 256G G 2 /45), 

a 6 = 8(-2G /35 + 64Gg/27 - 3584G§/81 + (29/35 - 
128Gg/3)G 2 /3-4G 4 /21)/3 

a 7 = V2(-5/896 + 769Gg/1400 - 286G^/15 + 

77792G[j/243 + 32(-193/35 + 176G§)G G 2 /45 + 
272G^/75 + 736G G 4 /315) 



?(o 



c(2 



c(3 



c(4 



?(5 



Coefficients of £q in (|6£|) 

= 1 

= -2V2Go 

= 3/20 + 28G§/3 

= -2V2(Go/5 + 320Gg/27 + 8G 2 /15) 

= -3/224 + 9Gg/5 + 1144G^/9 + 64G G 2 /5 

= 2V2(G /28 - 16<3g/9 - 14336Gg/81 + 
4(-l/35 - 64Gg/9)G 2 - 4G 4 /35) 



K = 1/2 

Ki = -2V2Go 

K 2 = 44G2/3 

K 3 = (G - 944G3/3 - 16G 2 )/3\/2 

K 4 = 2(-7Gg/5 + 5008Go=/27 + 304G G 2 /15) 

K 5 = -3G /80 + 734Gg/45 - 35216Gg/27 + 2(1 - 
5008Gg/15)G 2 /3 - 16G 4 /15 

K 6 = 8(llGg/35 - 8236G^/135 + 277600Gg/81 + 
4(- 151/7 + 27016G§/9)G G 2 /15 + 52G^/5 + 
712G G 4 /105)/3 

Coefficients 9 t of 9 in (0) 

80 = 1/8 

01 = 3\/2Go/2 

2 = 5G2 

9 3 = V2(-G + 16(-7G§ + G 2 ))/8 

6 A = (-51Go + 16(535G|] + 69G 2 ))G /90 

6 5 = V2(81Go + 32(355Gg - 9(5G 2 - 8G 4 ) - 

40(743G ! + 130G 2 )Gg))/2880 

6 6 = 2(39Gg + 84(43G^ - 227G^) + 32(33985Gg + 

8001G 2 )Gg - 12(107G 2 - 138G 4 )G )/945 
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Coefficients 0% of 9 in ([roj), for fermions without spin 

O = e 1 = o 

02 = 16G§ 

6> 3 = -160\/2Gg/3 

4 - 8(-G§ + 140G^ + 16G G 2 )/3 

5 = 8V2(91Gg - 7360Gj^ - 1264G^G 2 )/45 

6 = 32(Gg/5 - 137G^/3 + 23176Gg/9 + 4(-l + 

452G§/3)G G 2 + 8G| + 24G G 4 /5)/9 

7 = 2V2(-613Gg/70 + 10048G^/9 - 431744G£/9 + 

16(1103/35 - 8032Gg/3)GgG 2 /3 - 1408G G|/3 - 
5504GgG 4 /35)/3 

Our results can be generalized for the case of spinless 
fermions. Here instead of Eq. ( p3|) and Eq. ( |35| ) one has 
to use 

K = ir X v F , 9=UVK--^=\ • (71) 
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